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ABSTRACT 

The mass domain where massive extrasolar planets and brown dwarfs overlap is still 
poorly understood due to the paucity of brown dwarfs orbiting close to solar-type stars, the 
so-called brown dwarf desert. In this paper we collect all of available data about close brown 
dwarfs around solar type stars and their host stars from literature and study the demographics 
of the brown dwarf desert. The data clearly show a short period and a medium mass gap in 
the brown dwarf period-mass distribution diagram (35 < msini < 55Afj up and P < 100 
days), representing the "driest land" in the brown dwarf desert. Observation biases are highly 
unlikely to cause this gap due to its short period and medium mass, of which brown dwarfs 
can be easily detected by previous RV surveys. Brown dwarfs above and below this gap have 
significantly different eccentricity distribution, which not only confirms that this gap is real, 
but also implies that they may have different origins. Our further statistical study indicates 
that brown dwarfs below this gap may primarily form in the protoplanetary disk through disk 
gravitational instability, while brown dwarfs above this gap may dominantly form like a stellar 
binary through molecular cloud fragmentation. Our discoveries have offered important insights 
about brown dwarf formation mechanisms and their possible relationships with planet and star 
formation. 
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1 INTRODUCTION 

Brown dwarfs (BD) are in the mass range of approximately 
13 — 80 Jupiter masses, having sufficient masses to burn deu- 
terium but not enou g h to burn hydrogen in th e ir inner cores 
jBurrows et all Il997l; Ichabrier & Baraffel l2000t burrows et ail 
2001; Spieg el. Burrows. & M ilsom 2011). The first discovery of a 
bona-fide BD (Rebolo et al. 1995; Nakajima et al. 1995; Oppen- 
heimer et al. 1995; Basri et al. 1996; Rebolo et al. 1996) came in 
the same year as the disc overy of the first extra-so lar planet around a 
solar type star, 51 Peg b I IMavor & Oueloall995h . One of the major 
achievements of high-precision radial velocity (RV) surveys over 
the past two decades is the identification of a brown dwarf desert, a 
paucity of brown dwarf companions relative to planets within 3 AU 
around main-sequence FGKM stars (Marcy & Butler 2000; Grether 
& Lineweaver 2006). Although the induced reflex RV signal by a 
close BD companion on a solar type star is well within the detec- 
tion sensitivities of the high precision RV surveys (~ 3 — 10 m/s), 
only a few dozens are known (Reid & Metchev 2008; Sahlmann 
et al. 2011a and references therein) compared to over five hundred 
giant planets detected so far by RV technique. The California & 
Carnegie Planet Search measured an occurrence rate of 0.7% ± 



0.2% from their sample of ~ 1000 target stars llVoet et al.ll2002l : 
|Pateletalj|2007l) . and the McDonald Observatory Planet Search 
shows a similar rate of 0.8% ± 0.6% from a search sample of 250 
stars dWittenmver et al.l 120091) . Sahlmann et al. (2011a) obtained 
an upper limit of 0.6% for the frequency of close BD companions 
based on the uniform stellar sample of the CORALIE planet search, 
which contains 1600 solar type stars within 50 pc. 

To assess the reality of the brown dwarf desert, 
iGrether & Lineweaverl d2006l) performed a detailed investiga- 
tion of the companions around nearby Sun-like stars. They found 
that approximately 16% of nearby Sun-like stars have close (P < 5 
yr) companions more massive than Jupiter: 11% ± 3% are stellar 
companions, < 1% are BDs, and 5% ± 2% are gi ant planets. 
Althou gh the close BDs are rare around solar type stars. iGizis et alj 

as rare at wide separations 
_ as at close separations. 
Lafreniere et al. d2007h obtained a 95% confidence interval of 
1.9~t\ g% for the frequency of 13— 75Mj up companions between 
25— 250AU in the Gemini Deep Planet Su rvey around 85 nearby 
young stars. iMetchev & Hillenbrand (2009) inferred the frequency 
of BDs in 28—1590 AU orbits around young solar analogs is 
3.2^2 ^% from an adaptive optics survey for substellar companions 
around 266 Sun-like stars. 



rtitnougii uie close dus aie rare aiounu so 
d200ll) su ggests that BDs might not be a s 
(see also [M etchev & Hillenbrand 2004) 
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BDs are traditionally believed to form like stars, through 
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gravitationa l collapse and/o r fragm e ntation of molecular 
cloud s JPadoan & Nordlundl 12004 iHennebelle & Chabried 
2008). A recently found self-gravitating clump of gas and 
dust has a mass (0.015 — 0.03 Mp,) in the BD regime 
jAndre. Ward-Thompson. & Greavesl |2012h . which supports 
the idea that BDs could form like a sta r. On the other ha nd, 
companions with masses up to 10 Mj up ( Alibert et ah! 120050 or 
even 38 Mj up (Mordasini et al. 2009) may form in protoplanetary 
disks according to the standard core-accretion planet formation 
theory. Because of this, the brown dwarf desert is commonly 
interpreted as the gap between the largest mass objects that can 
be formed in protoplanetary disks, and the smallest mass clumps 
that can collapse and/or fragment in the vicinity of a protostar. 
The mass function of close stellar companions shows a linear 
decrease in log(M) toward the BD mass range from both stellar 
mass and planetary mass directions (Grether & Lineweaver 2006). 
In comparison, the mass function of isolated substellar objects 
seems to be roughly flat in log(M) down to masses ~ 20 Mj up , 
both in the field and in clusters (Luhman et al. 2000; Chabrier 
2002). This indicates that close BD companions may form in a 
different way from those formed in the field and clusters. 

As such, statistical properties of close BD companions as well 
as how these properties are related to their host stars, contain a lot 
of information about the poorly understood BD formation mech- 
anisms and their relationships with star and planet formations in 
close orbital environments. These statistics may also be important 
to investigating how additional important processes such as tidal 
evolut ion and disk-planet intera c tion affect close BD pr operties 
(e.g.. lArmitaee & Bonnelil 120021; iMatzner & Levir] l2005t) . Given 
that the close BD occurrence rate is < 1%, currently there has yet 
been a single large, relatively uniform RV survey capable of produc- 
ing a large homogeneous sample of BD companions for a meaning- 
ful statistical study (Marcy et al. 2000; Ge et al. 2008; Sahlmann et 
al. 201 la). However, all of the previous RV planet surveys have suf- 
ficient RV sensitivity and time baseline (at least more than 2 years, 
e.g., Ge et al. 2008; Ge et al. 2009; Eisenstein et al. 2011) to de- 
tect a majority of close BDs around solar type stars due to the large 
RV amplitude (on the order of ~ 1000 m/s, vs. a few to a few tens 
m/s RV precisions in previous RV planet surveys). It is, therefore, 
possible to combine close BD companion samples together for a 
statistical study without major biases. 

In this paper we assembled a catalog of all the BD companions 
discovered around solar type star from literature and used it to con- 
duct a statistical study. We have found tentative evidence for the ex- 
istence of two different populations of BD companions. We present 
the BD catalog assembled in this study in §2 , statistical properties 
of BD companions in §3 and discuss their implication for different 
BD formation scenarios in §4. We summarize our main results in 
§5. 



2 CATALOG DESCRIPTION 

We have collected data from literature about the currently known 
BD (candidates) companions around FGK type stars. Most of them 
have known Kepler ian orbits, except for H IP 78530 (discovered 
by direct imagin g; Lafren iere et al.l 1201 ll) . KOI-205.01 (unpub- 
lished vet; | Santerne et al.ll2012h and KOI-554.01 (unpublished yet; 
ISanterne et al.l2012h . Properties of the BDs (minimum mass, pe- 
riod and eccentricity) and their host stars (mass, effective temper- 
ature, surface gravity and metallicity) are summarized in Table [T] 



log P (days) 

Figure 1. Period distribution of known BD companions around solar type 
stars. 




io' icr 

Period (d) 

Figure 2. Cumulative mass distribution of brown dwarf candidates. Three 
lines with three RV precisions, 50 m/s, 100 m/s and 150 m/s, are also shown. 



For those which are transiting their parent stars or have astrometry 
measurements, true masses are also given besides Msini. 



3 OBSERVED PROPERTIES OF BROWN DWARFS 
3.1 Orbital Period Distribution 

The distribution of orbital periods of BDs has two main features 
(Fig.Q3 : a relatively flat distribution inside P ~ 100 days and a 
sharp jump beyond P ~ 100 days. The drop beyond P ~ 1000 
days is likely due to the observational incompleteness since 1) 
it is more difficult to detect a BD companion over a long pe- 
riod than a short period with RVs and 2) some RV surveys (such 
as the SDSS-III MARVELS, Ge et al. 2008; Ge et al. 2009; 
Eisenstein et al. 2011) do not cover beyond this period. It is ev- 
ident that the number of BDs increases with the orbital period, 
even though RV and transit observations are biased toward dis- 
covering objects in short-periods. The position of the maximum 
of the distribution is unknown due to the different duration limit 
of most of the old surveys (several thousand days). This increas- 
ing distribution is consistent with the results from high contrast 
and high angular resolution imaging surveys, which find evidence 
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Figure 3. Cumulative mass distribution of brown dwarf companions. Brown 
dwarfs with periods greater and less than 100 days are shown as dashed and 
solid lines respectively. 



for a higher fraction of BP companions at wide orbits than at 
close orbits l Lafreniere et all 2007: Metc hev & Hillenbrand 120091 : 



IChauvin et al]l20ld : ljanson et al.ll2012h 



JCumrnine. Marcv. & Butlei 


1999; Udrv. Mayor, & Santos 


120031: IMarcv et alJ 120051: 


Udrv & Santos 


20071) and binaries 


jDuauennov & Mayor! Il991 


) show an increasing number distri- 



bution with period. However, there is no 3-day pile up with the 
BD distribution as shown in giant planets (Udry et al. 2003). The 
reason may be that BDs initially form in the protoplanetary disks 
further away from the host stars where protoplanetary disks have 
more materials to efficiently form BDs than at the close-in regions, 
and the migration mechanism may not efficiently move such a 
massive body to a short period orbit (Trilling et al. 1998; Nelson 
et al. 2000; Trilling, Lunine & Benz 2002). On the other hand, 
BDs forming at the same time as the primary stars may migrate 
inward quickly in the initi al gas rich disks and thus be destroyed 
via mergers with the stars (Armitage & Bonnell 2 0021) . 



3.2 The Period-Mass Diagram 

The period-mass distribution of BDs shows a statistically signif- 
icant gap at the short period and medium mass region as illus- 
trated in Fig. |2 In this plot, a rectangular area with P < 100 d 
and 35M j up < M < 55Mj up is highlighted to show the gap 
within which BDs are nearly depleted while there are numerous 
BDs around this region. It appears that this gap is real since it is 
unlikely caused by detection sensitivity (RV precision) or observa- 
tion biases (due to survey incompleteness). Previous RV observa- 
tions have detected many BDs with masses less than the lower limit 
(~35 Mj up ) of this region and also BDs with periods significantly 
longer than the period limit of this region (~ 100 days). RV sensitiv- 
ities with three moderate RV precisions, 50 m/s, 100 m/s and 150 
m/s, shown in Fig. [2] clearly illustrate that any of BDs in the gap 
should be detected easily with these moderate RV precisions. To 
further verify if this feature is real, we divided the BD sample into 
two groups according to their periods (P < 100 days and P > 100 
days) and plotted their mass cumulative histograms in Fig. [3] for 
comparison. It is apparent that a depletion of BDs with masses be- 
tween 35 and 55 A/j up appears in the cumulative histogram for the 
short period group while no depletion of BDs appears in the cumu- 



lative histogram for the long period group. We carried out a sim- 
ple Monte Carlo experiment to test the emptiness of this gap on 
the period-mass diagram. There are a total of 25 BDs with period 
shorter than 100 days in our BD sample. We assumed simply that 
their masses are uniformly distributed between 13 and 80Mj up . 
Then we drew their masses randomly from this uniform distribu- 
tion and counted how many of them will fall in the mass range of 
35 to 55Mj up . We found the probability that less than 3 BDs will 
fall in this gap is 0.9%, which corresponds to a 2.6a significance. 
A larger BD sample in the future will be better to assess the signif- 
icance of this gap. 

The appearance of this depleted region in the BD period-mass 
diagram has naturally divided BDs into two mass groups: one with 
masses greater than 42.5Mj up and the other less than 42.5Mj up . 
Their properties and origins may be different. We further explore 
properties of these two groups and study possible origins. 

3.3 Orbital Eccentricity Distribution 

The orbital eccentricities show great difference for the two BD 
groups with masses greater and lower than 42.5A/j up , respectively. 
Fig.f4]shows the period-eccentricity distribution of all known BDs. 
The period-eccentricity distribution of BDs with masses greater 
than 42.5M.7up is consistent with a circularization limit of ~ 12 
days, which is simila r to that found in nearby stellar binaries 
lRagh avanet"aT]|2010l) . It is clear that there are a significant num- 
ber of BDs with 300 d < P < 3000 d and e < 0.4 for BDs 
with masses lower than 42.5A/j up , but no BDs with masses greater 
than 42.5A/j up . We have done a two-dimensional Kolmogorov- 
Smirnov (K-S) test for the period-eccentricity distribution of BDs 
with masses greater and lower than 42.5Mj up . The probability that 
these two BD samples are drawn from the same distribution in the 
period-eccentricity plane is 1.7%. 

Next we are going to compare the per iod-eccentricity d i stribu - 
tion of the BDs to that of stellar binaries. Halb wachs et al.1 d2003l) 
have studied the statistical properties of a sample of 89 FGK type 
main-sequence binaries with periods up to 10 year. Here we choose 
to use their 89 binary sample for our comparison. To compare the 
period-eccentricity distribution between BDs and stellar binaries 
we made use of a two-dimensional K-S test. The probability for the 
period-eccentricity distribution to be the same is 18% between BDs 
with masses above 42.5Mj up and the stellar binary sample, and 

0. 1% between BDs with masses below 42.5Mj up and the stellar 
binary sample. These results suggest that BDs with masses greater 
than 42.5Mj up have a very similar period-eccentricity distribution 
to that of stellar binaries, and their formation mechanisms may be 
similar. 

The difference of the orbital eccentricities for BDs with differ- 
ent masses is further illustrated in the mass-eccentricity plot shown 
in Fig. [5] By including all currently known planets (from exo- 
planet.org) and brown dwarfs (this paper) in this plot, a clear trend 
is shown: all the known giant planets and BDs with masses below 
~ 42.5Mj up have the eccentricity distribution following a trend, 

1. e., the more massive the giant planet/BD is, the lower maximum 
eccentricity it tends to have while BDs above this mass threshold 
do not show such a trend, instead showing more diversity in their 
eccentricities. 

3.4 Metallicity of the BD Host Stars 

The BD host stars in this study have a mean metallicity of [Fe/H] 
= —0.04 with a standard deviation of 0.28. For comparison, 
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Figure 4. Period-eccentricity distribution of brown dwarf candidates. Brown 
dwarf candidates with masses above and below 42.5Afj up are shown as 
circles and diamonds. Brown dwarf candidates with true masses measured 
using transiting observations or astrometry measurements are shown as 
crosses. 
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Figure 5. Mass-eccentricity distribution of exoplanets and brown dwarf 
candidates. The solid line shows a trend that the more massive the gi- 
ant planet/BD is, the smaller maximum eccentricity it tends to have, 
which breaks at ~ 42.5Afj up . The vertical dotted line shows Msini = 
42.5Afj up . Planets are shown as squares. Brown dwarfs with masses above 
and below 42.5Mj up are shown as triangles and circles respectively. 



Raghavan et al. I d2010h has carefully studied a sample of 454 nearby 
solar type stars, which have a mean metallicity of —0.14 with a 
standard deviation of 0.25. Six stars in that sample have BD com- 
panions, which have a mean metallicity of [Fe/H]= —0.05. The 
mean metallicity of our BD sample is slightly higher than the mean 
metalli cally of volume limited nearby F GK dwarf stars. For ex- 
ample, [FavaTaj^ic^la^^SclOTtirio] |l993) have analyzed a volume 
limited sample of 91 G and K dwarfs, yielding a m ean metallicity of 
lFe/Hj = —0.08 with a standard deviation of 0.26. Nordstro m et al.l 
(2004) have derived metallicity for 16682 nearby F and G dwarf 
stars with a mean of —0.14 and a dispersion of 0.19 dex. Sousa et 
al. (2011) found that a mean metallicity for the CORALIE survey 
sample of 1248 stars and the HARPS survey sample of 582 stars is 
[Fe/H]=— 0.11 and —0.10, respectively. However, since some exo- 
planet surveys choose samples biased towards metal rich stars (e.g., 
Valenti & Fischer 2005), the slightly higher mean metallicity of BD 
host stars is possibly caused by the sample bias. 

After comparing the BD host star met allicities with the vol- 
ume limited sample from lSousa et al.l d201 lh and the planet search 
sample from Valenti & Fischer (2005), we find that the BD host 
star metallicity distribution is consistent with the combination of 
these two samples. This means we cannot interpret the BD host star 
sample as metal rich. We compared metallicities of BD host stars 
with that of giant planet (lMj up < msini < 5M,j up ) host stars. 
The data of g iant planet host star s is taken from the Exoplanet Or- 
bit Database JWrightet aT1l201ll) . A K-S test shows that the prob- 
ability of these two samples selected from the same distribution is 
2 x 10~ 4 . This indicates that the two samples are significantly dif- 
ferent from each other. 

We investigated the correlation between BD host star metallic- 
ities and BD masses. The Spearman's rank correlation coefficient 
of the BD mass and their host star metallicity is 0.07 with a 61% 
significance, suggesting there is no significant correlation between 
the BD mass and their host star metallicity. 

We also divided the BD companions into two subsample ac- 
cording to their masses. The cumulative metallicity distribution for 
host stars with BD companion masses above and below 42.5Mj up 
is shown in Fig. [6] The main difference between these two distribu- 
tion is at the lower metallicity end. Currently no BDs with masses 
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Figure 6. Cumulative metallicity distribution of host stars with BD com- 
panion masses above (dashed line) and below (solid line) 42.5Afj up . Also 
shown for comparis on are cumulative metal licity distribution of the planet 
search s ample fromlValenti & F ischer ( 2005 ), the CORALIE planet search 
sample ( Sousa et al. 201 1 ) and giant planet sample from the Exoplanet Orbit 
Database iwright et alj201 lh . 

greater than 42.5Afj up have been found around stars with [Fe/H]< 
—0.5, while several BDs with masses lower than 42.5Afj up have 
been found in this metallicity regime. A K-S test has been con- 
ducted to show that the probability of these two sample metallicities 
are drawn from the same distribution is 70%, which suggests there 
is no significant difference between these two samples regarding 
their metallicity distribution. 



4 DISCUSSION 

4.1 Two Different Brown Dwarf Populations 

Our study suggests that BDs with masses lower than ~ 43Mj up 
have an eccentricity distribution consistent with that of giant plan- 
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ets in the mass-eccentricity diagram while BDs with masses above 
~ 43Mj U p have the star-like eccentricity distribution (Figure |3]. 
Our mass limit is consistent with the minimum of the mass func- 
tions of planet and stel lar companions in the BP m ass region, 
43I23MJUP, derived by iGrether & Lineweaverl d2006l) using their 
stellar companion and giant planet sample within 50 pc around 
the sun. This mass function minimum is also consistent with 
that derived by ISahlmann et al.l d2011bl) . who found a void in the 
mass range between 25 and 45 Mj up using the data from the 
CORALIE radial-velocity survey. They further suggested that there 
may be a possible d ividing line between m assive planets and sub- 
stellar companions. ISchneider et al.l (1201 lh has chosen arbitrarily 
and probably provisionally 25M,t up as t he upper limit of mas- 
sive planets based on previous stud ies (e.g.. lSahlmann et alj2011bl 
iBaraffe. Chabrier. & BarmanfcoiOl) . 

Our BD sample is, therefore, naturally divided into two dif- 
ferent groups with the mass limit of ~ 42.5M.j up . The eccentric- 
ity distribution of low mass BDs appears to be consist ent with the 
prediction from the "planet-planet scattering" mo del dRasio & Ford 
ll996tlFord & Rasidl2008l ; IChatteriee et al.ll2008l) . while the eccen- 
tricity distributio n of massive BDs app ears to be similar to that of 
stellar binaries (Halbwachs et al. 2003) as shown in ^ 13.31 The ec- 
centricity distributions of these two groups support that BDs may 
form differently: BDs below this mass limit form in protoplanetary 
disks around host stars, and above this mass limit form like stellar 
binary systems. This is supported by our analysis results. The ex- 
istence of a large population of long period low eccentricity BDs 
(P > 300 days and e < 0.4) serves as evidence to support the BD 
formation scenario in the protoplanetary disks for those compan- 
ions with masses below 42.5Afj up . While the lack of long period 
and low eccentric BD companions with masses above 42.5Afj up 
appears to support the BD formation scenario like stellar binary 
formation. Nevertheless, a small number of BDs in each of these 
two mass regions may form in an opposite formation mechanism, 
but our sample is not sufficient to distinguish these minor groups. 



4.2 BD formation mechanisms 

Our study of metallicity distrib ution of B Ds appears to support 
the disk instability mechanism dBossll 19971) for those BDs formed 
in the planetary disks, while it is inconsiste nt with the predic- 
tion of core-accretion formation mechanism dPollack et all 1 19961 : 
llda & Linll2004l : lAlibert et al .1120051) . Currently, there are two hotly 
contested theories about giant planet formations: core-accretion 
and disk gravitational instability. In the core-accretion scenario, 
giant planets form more efficiently around metal rich planetary 
disks than metal poor planetary disks because the higher the grain 
content of the disk is, the easier the metal core i s to build for 
giant planets dPollack et al] 1 19961 : llda & LirJ |2004 lAlibert et alj 
120051) . In contrast, the disk gravitational instability process al- 
lows similar formation efficiency for both metal rich and metal 
poor giant planets. Previous results show strong correlation be- 
tween metalli city and giant planet occurrenc e rate around so- 
lar type sta rs (Santos, Israe lian. & Mavonl200ll: IValenti & Fischen 
l2005l : IJohnson et alj |20 1 0) strongly support th e core-accretion for- 
mation scenario for giant planet f ormation dPollack et al] 1 1 99*61 : 
llda & Linll200l : lAlibert et al.ll2005l) . However, our result, showing 
no correlation between the BD occurrence rate and metallicity, ap- 
pears to support the disk instability mechanism for formation of 
BDs in planetary disks. 

For stellar companions, there appears to be a weak anti- 
correlation between metallicity and stellar companion occurrence 



rate (Ragh avanetail hoiO). i.e., lower-metallicity clouds might be 
more likely to fragment to form binary stars. If this applies to mas- 
sive BDs, then we expect this anti-correlation. However, our current 
sample is too small to tell, but is at least not inconsistent with this 
statistics. 



4.3 Upper Mass limit of Planets 

Our study appears to help to address the upper mass limit of "plan- 
ets". Here we define a "planet" as an object formed in a protoplan- 
etary disk from the companion formation point of view. It has been 
long known that the low mass end of brown dwarfs overlaps the 
high mass end of massive planets formed in protoplanetary disks. 
From theoretical si mulations, planets with masses up to 10 Mj up 
dAlibert et ai]|2005l) or even 38 Af Jup (Mordasini et al. 2009) may 
form in protoplanetary disks. However, it is difficult to determine 
the upper mass limit of planet from observations since the high 
mass planets detected may actually be brown dwarfs formed as low- 
mass binaries. One way to tackle this problem is to search for de- 
bris disk around massive planet/BD host star, which support the 
idea th at these objects are formed like a planet. iMoro-Martm et al. 
d20ld) have found evidence of debris disk around two massive 
planet host star (HD 38529 and HD 202206), which supports for- 
mation of planets up to 17 Mj up in a protoplanetary disk. The other 
way to resolve this problem is through the study of mass spectrum 
of massive planets/BD. Motivated by the observed mass distribu- 
tion, this planet upper mass l imit has been se t at around 43l,t Afj up 
(Grether & Lineweaver 2006) or 25 Mj up (iSah lmann et al]|201 lbl : 
ISchneider et al]|201 ll) . Analysis conducted in this paper was also 
able to provide a clue about this upper mass limit. The cumula- 
tive BD mass distribution (Fig.|3j suggests this limit is in the range 
30 — 60 Af.j up and mass-eccentricity relation (Fig. [5} suggests 
that this limit is around 43 Mj up . It appears that the upper limit 
is around 43 A/j up , which i s consistent with the result found by 
iGrether^&I^meweavei d2006l) , but is slightly larger than that from 
i Sahlmannetal .11201 lbh . 



5 SUMMARY 

We have searched literature and presented a catalog of BD compan- 
ions around solar type stars found by radial velocity, transiting and 
astrometry observations. We have studied distribution of different 
parameters of BD companions around solar type star, and found 
that: 

(1) BD companions have an increasing distribution with pe- 
riod, similar to giant planets and low-mass binaries; 

(2) BD companions are almost depleted at P < 100 d and 
30Mj up < M < 55Mj up in the period-mass diagram; 

(3) BD companions with masses below 42.5Mj up have eccentricity 
distribution consistent with that of massive planets; 

(4) BD companions with masses above 42.5A/j up have eccen- 
tricity distribution consistent with that of binaries, which shows 
the expected circularization for periods below 12 days, caused by 
tidal forces over the age of the Galaxy, followed by a roughly flat 
distribution; 

(5) Host stars of BD companions are not metal rich, and have 
significantly different metallicity distribution when comparing with 
host stars of giant planets, suggesting a formation scenario at least 
partly different from the core-accretion scenario. 

The distribution of BD and their host star properties presented 
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in this paper may lend support to such a picture: (1) BD compan- 
ions with masses below 42.5Mj up form in a protoplanetary disk 
through the the instability-fragmentation scenario, and their eccen- 
tricity is excited through scattering with other objects formed in this 
disk or interactions with disk/third body; (2) BD companions with 
masses above 42.5Mj up form like stars, through molecular cloud 
fragmentation, similar to the formation of a stellar binary system. 
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-0.27 ± 0.02 


48,53 


' 1 Bryant 



Table 1 (cont'd) 



Object 


Mo 


Aie sin i 


Period 


Eccentricity 




T et f 


log(s) 


[Fe/H] 


References 




(M Jup ) 


(M Jup ) 


(day) 




(M ) 


(K) 


(cgs) 








DO. f ± J.U 




26772 _ 1059.2 




1 .Uo ± U.U4 




4.42±0.06 


U.00 ± U.O / 




nly I 1 JO / J 




O / . O _J_ jS. a 


i jdd.o _c o . ( 


U. J / o Hz U.UUQ 


2 7 ± 3 




z. jy zn u . x u 


n 1 4 4- n in 


SO 


HD 136118 


42 Il8 


12.0±0.47 


1 188.0±2.0 


0.34±0.01 


1.24±0.07 


6097 ±44 


4.16±0.09 


-0.01±0.05 


59, 60 


HD 217786 


13.0±0.8 


1319 ±4 


0.40±0.05 


1.02±0.03 


5966±65 


4.35±0.11 


-0.135±0.043 


61 


WASP-30 


60.96 ± 0.89 




4.156736±0. 000013 


(adopted) 


1.17 ± 0.03 


6201 ±97 


4.28±0.01 


-0.03 ± 0.10 


13 


COROT-15 


63.3 ± 4.1 




3.06036±0. 00003 


(adopted) 


1.32 ± 0.12 


6350±200 


4.3±0.2 


0.10 ± 0.20 


14 


LHS 6343 


62.9 ± 2.3 




12.71382±0. 00004 


0.056± 0.032 


0.370 ± 0.009 




4.851±0.008 


0.04 ± 0.08 


15 


Corot-3b 


21.66 ± 1.00 




4.25680±0. 000005 


(adopted) 


1.37 ± 0.09 


6740±140 


4.22±0.()7 


-0.02 ± 0.06 


19 


X0-3b 


11.8 ± 0.59 




3.1915239±0. 0000068 


0.26±0.017 


1.213 ± 0.066 


6429±100 


4.244±0.041 


-0.177 ± 0.080 


20 


KOI-423b 


1M +o.» 




2 1.0874 ±0.0002 


0.121+ 022 


, |+0.07 


6260±140 


4.1±0.2 


-0.29 ± 0.10 


23 










-0.023 


-0.06 




4 229+ ' 012 
-0.019 






KELT- lb 


27 23~*" " 5 
z/ z -0.48 




1.2175 14±0. 000015 


0099 +0 " 01 

0.0069 


1.324 ± 0.026 


6518±50 


0.008 ± 0.073 


56 


KOI-205.01 


35 




11.72 






5060 


4.57 


-0.17 


57 


KOI-554.01 


80 




3.66 






5835 


4.64 


-0.08 


57 
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